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The chiral phase transition is investigated within the framework of thermal field theory using 
the linear sigma model as an effective theory. We concentrate on the meson sector of the model, 
and calculate the thermal effective potential in the Hartree approximation by using the Cornwall- 
Jackiw-Tomboulis formalism of composite operators. The thermal effective potential is calculated 
for A'^ = 4 involving as usual the sigma and the three pious, and in the large- A'^ approximation 
involving A*' — 1 pion fields. In the iV = 4 case, we have examined the theory both in the chiral 
limit and with the presence of a symmetry breaking term, which is responsible for the generation of 
the pion masses. In both cases, the system of the resulting gap equations for the thermal effective 
masses of the particles has been solved numerically, and we have investigated the evolution of the 
effective potential. In the A'^ = 4 case, there is indication of a first-order phase transition, and the 
Goldstone theorem is not satisfied. The situation is different in the general case where we have 
used the large-A'^ approximation. The Goldstone theorem is satisfied, and the phase transition 
appears as second-order. In our analysis, we have ignored quantum fiuctuations and have used the 
imaginary time formalism for calculations. We extended our calculation in order to include the full 
effect of two loops in the calculation of the effective potential. In this case, the effective masses are 
momentum dependent. In order to perform the calculations, we found the real time formalism to 
be convenient. We have calculated the effective masses of pions at the low-temperature phase and 
we found a quadratic dependence on temperature, in contrast to the Hartree case, where the mass 
is proportional to temperature. The sigma mass was investigated in the presence of massive pions, 
and we found a small deviation compared to the Hartree case. In all cases, the system approaches 
the behaviour of the ideal gas at the high temperature limit. 
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is the order of this transition, is it first-order with latent 
heat, or second-order, or maybe a crossover between the 
phases? Lattice calculations suggest that when we con- 
sider two massless quarks, the transition is second order 
and the same suggest other approaches based on effective 
models. If we consider three massless flavours of quarks, 
the transition is probably first-order IB| . The aim of 
this work is to study the chiral phase transition. 
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I. CHIRAL SYMMETRY BREAKING AND 
RESTORATION 

A. Overview 

The study of matter at very high temperatures and 
densities and of the phase transitions which take place 
between the different phases, has several very interesting 
aspects. It has been the subject of intense study dur- 
ing the last few years, because of its relevance to particle 
physics, astrophysics and cosmology. According to the 
standard big bang model, it is believed that a series of 
phase transitions happened at the early stages of the evo- 
lution of the universe, the QCD phase transition being 
one of them 13, 0, S H S • 

There is hope at present that it may also be possible 
to probe this transition in the laboratory in experiments 
involving relativistic heavy ion collisions. Experiments 
of this type carried out at CERN may have reached the 
phase transition 0. Further experiments are planned in 
the future, and the results could possibly improve our 
knowledge on subjects such as the restoration of chiral 
symmetry, the nature of the quark gluon plasma and the 
color superconducting phase, as well as the physics of 
neutron and quark stars Q . 

There are two main issues related to the QCD phase 
transition, namely restoration of chiral symmetry (chi- 
ral phase transition), and deconfinement of quarks and 
gluons to form the so called quark-gluon plasma. At 
present, it is not at all clear what the relation is between 
these two phase transitions, if they happen at the same 
temperature, or if they are independent. Another issue 



proximate Goldstone bosons. It is generally believed that 
at some temperature, or baryon density, the chiral sym- 
metry could be restored. 

An important part when studying questions like the 
restoration of spontaneously broken symmetries is the 
construction of order parameters, which characterise the 
way in which the symmetry of the system under con- 
sideration is realised. These quantities are zero in the 
phase where the symmetry is manifest but non zero in 
the spontaneously broken phase. A classic example of an 
order parameter is the magnetisation of a ferromagnetic 
substance, which is non-zero below the Curie tempera- 
ture, but disappears at higher temperatures. The sys- 
tem undergoes a transition from an asymmetric ordered 
state, with non-zero magnetisation at low temperature 
to a symmetric disordered state, with zero magnetisation 
at temperatures well above the Curie point. We usually 
encounter two types of phase transition. In first-order 
transitions, the order parameter jumps discontinuously 
from its value in one phase to that in the other (usually 
zero). In contrast, during second order transitions, the 
order parameter vanishes continuously 0, |^ . 

An important order parameter for the chiral phase 
transition of QCD is the quark condensate, a measure 
of the density of quark-antiquark pairs that have con- 
densed into the same quantum mechanical state. They 
fill the lowest energy state - the vacuum of QCD - and 
as a result the chiral symmetry is broken, since there is 
no invariance under chiral transformations [lol IT^ . !l3| . 
It is expected that at very high temperatures, the quark 
condensate disappears, and the system is chirally sym- 
metric. 

There are two main paths to study the chiral phase 
transition, namely lattice QCD methods and effective 
field theories. Lattice QCD can tell us many things 
about the phase transition, but cannot be used to study 
its dynamics. Hence, the need for models which do not 
suffer this restriction. In the case of chiral symmetry, 
a model with the correct chiral properties is the linear 
sigma model, a theory of fermions (quarks or nucleons). 
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interacting with mesons This model has been used 
extensively as an effective theory in the low energy phe- 
nomenology of QCD, describing the physics of mesons, 
and it is well suited for a study of the chiral phase tran- 
sition [lol |. In particular the model is very popular 
in studies of the so-called disoriented chiral condensates 
(DCC). We review the model and discuss how it is related 
to chiral symmetry in Section llBI For this presentation 
we have followed closely the lecture notes on spontaneous 
breaking of chiral symmetry by Li 

The appropriate framework for study of phase transi- 
tions is thermal field theory or finite temperature field 
theory, a combination of quantum field theory and sta- 
tistical mechanics. Within this framework, the finite- 
temperature effective potential is an important and well- 
used theoretical tool. The use of such techniques goes 
back to the 1970's when Kirzhnits and Linde [iM 
first proposed that symmetries broken at zero tempera- 
ture could be restored at finite temperatures. Subsequent 
work by Weinberg 17], Dolan and Jackiw as well as 
many others, resulted in a wide adoption of the effective 
potential as the basic tool in such studies. The basic 
ideas about the notion of the effective potential are re- 
viewed in the Section ITCl However, a detailed analysis 
can be found in the original pap ers, and also in recent 
textbooks on the subject [TllM l2lll23| . 

The finite-temperature effective potential V{(j), T) is 
defined through an effective action r(0), which is the gen- 
erating functional of the one-particle irreducible graphs, 
and it is related to the free energy of the system [18l |. 
A generalised version is the effective potential V(0, G) 
for composite operators introduced by Cornwall, Jackiw 
and Tomboulis (CJT) ^2^] and later formulated at finite 
temperature by Amelino-Camelia and Pi p3. l25l l26l| . 
The composite operator method in calculating thermo- 
dynamic potentials for many-body systems, was also in- 
troduced in the 1960s by Luttinger and Ward (see 
also 13), and has recently been used in studies of sys- 
tems in and out of equilibrium l30l| . 

When one is working in the context of field theory at 
finite temperature, there are basically three main paths 
to follow, namely the imaginary time formalism, the real 
time formalism and thermo-field dynamics p3 |. Which 
of these is the more convenient depends on the applica- 
tion at hand. There are advantages and disadvantages to 
each of the three methods. We do not use thermo-field 
dynamical methods in what follows, but we will use both 
real and imaginary time formalisms. We give a few basics 
about the two types of formalism in Section ITT?! Finally 
in Section ITeI we present the calculation of the effective 
potential for the \(tf^ theory, using the CJT method. 

The outline of the remaining sections is as follows: in 
Section |n] we apply the CJT method to calculate the ef- 
fective potential for the 0(4) linear sigma model in the 
Hartree approximation. We numerically solve the resul- 
tant system of gap equations, both in the chiral limit 
and with the presence of a linear term, which breaks the 
chiral symmetry of the Lagrangian. We have repeated 



these steps for the generalised version of the linear sigma 
model, with — 1 pion fields in the large-A^ approxi- 
mation, in Section IIIII . Then in Section IIVI we extend 
our Hartree approximation for the 0(4) model in order 
to include the effects of the so-called sunset diagrams. In 
these calculations, we find the real time formalism more 
convenient. 

In SectionEI we review recent attempts to incorporate 
the linear sigma model and/or the CJT formalism in or- 
der to investigate the physics of the hot pion gas, and the 
chiral phase transition. Finally, in Section lvTl we outline 
our conclusions. 



B. Chiral symmetry 

The symmetry principle is possibly one of the most im- 
portant ideas in the development of high energy physics. 
As it is well known, the symmetries of a physical sys- 
tem lead to conservation laws, which have as conse- 
quence many important relations for the physical pro- 
cesses. However, many of the symmetries which we ob- 
serve in nature are approximate symmetries rather than 
exact ones. A very interesting mechanism which is re- 
lated to these ideas is the spontaneous symmetry break- 
ing (SSB) and it seems that it has played an individual 
role in the development and our present understanding 
of particle physics. The main characteristic of sponta- 
neous symmetry breaking is the fact that it is related to 
the ground state of the theory. The idea of SSB was first 
appeared around 1960 in studies of superconductivity in 
solid state physics by Nambu and Goldstone 31,32]. One 
of the important consequences of SSB is the presence of 
massless excitations which are reffered to as the Nambu- 
Goldstone bosons. Later, these ideas were implemented 
in particle physics. In combination with SU{i) x SU{2>) 
current algebra, SSB has been quite successful in our 
understanding of the chiral symmetry, a important in- 
gredient of the low energy phenomenology of the strong 
interactions. However, although SSB has been quite suc- 
cessful in explaining many interesting phenomena, the 
implementation of this symmetry the theoretical frame- 
work has been done rather arbitrarily, and the origin of 
SSB is not absolutely clear. 

A popular model which implements the ideas of chiral 
symmetry is the linear sigma model which has a long and 
interesting history. It was originally constructed in the 
1960's, as a model to study the chiral symmetry in the 
pion-nucleon system 'l4]. Later the spontaneous sym- 
metry breaking and PCAC (partially conserved axial cur- 
rent) were incorporated. Although this model is not quite 
phenomenologically correct, it remains an attractive ex- 
ample with interesting features, which displays many im- 
portant aspects of broken symmetries. Nowdays it is a 
common belief that the strong interactions are the best 
described by QCD, however the linear sigma model, and 
also various quark models, still serve as effective theories 
at low energies, since at this regime due to confinenement 
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it is difficult to calculate directly from QCD. 

The Lagrangian for the linear sigma model is given by 



1 



(1) 



where the spinor field -0 is a massless isodoublet nucleon 
(or quark) field. The pseudoscalar field tt = (tti, 7r2, tts) 
is an isotriplet of pion fields and a is the isosinglet field. 

This Lagrangian is now clearly invariant under the (in- 
finitesimal) SU{2) transformations, 

Sijj — —ie ■ t/2tP, Stt — extt, 6a — , (2) 

and the associated current is 

J°^(a;) = (tt X 9^7r)° + ^V'T^r'^V' • (3) 

We can verify that the above Lagrangian is also invariant 
under a further set of transformations, namely 

6ip — —irj ■ t/2jqi/j, 6tt — rja, 5a = —ri-Tr. (4) 

Then we have another set of conserved currents 

J^f" = (aS^Tr"^ - 7r"9'^a) + ^ipY'lsT^i; , (5) 

and the corresponding charges are given by 

Qa^ j d?xr\x), Ql^ j cPxjf . (6) 



A point that should be made here, is that the transfor- 
mations given in Eq. Q), involve changes in parity and 
as a consequence the current in Eq. ||SJ) is an axial vec- 
tor and the associated charges are pseudoscalars. Using 
the canonical commutation relations, we can derive the 
algebra generated by those charges 



(7a) 
(7b) 
(7c) 



which show that the are the generators of an SU{2), 
and that the Qf's transform as an isovector under this 
SU{2). We can now consider the combinations Q1 and 
defined as 

Ql^liQ^ + Qt), QB = \iQ''-Qt) (8) 

to verify that they satisfy the following relations, 

[Ql,Ql]^teabcQl , (9) 



Therefore and commute with each other, and they 
generate separate SU{2) algebras. Then this combined 
algebra is called SUl{2) x SUr{2). 

There is also another way to describe the symmetry of 
the linear sigma model, that is the 0(4) symmetry which 
has the property of being isomorphic to SU{2) x SU{2). 
If we apply the transformation 



with 



(10) 



with R an orthogonal 4x4 matrix, then the scalar field 
= (cr, TTi, 772, TTs) trausforms as 4-dimensional vector, 

0» — * (f^i = Rij4>J — £ij(l>3 ■ (11) 

We then observe that the combination (pf = cr^ -I- tt^, is 
just the length of the vector and is clearly invariant 
under the rotations in 4-dimensions. If we take 



^3^ 



we find that 

7r' = 7r + Q;X7r + f3a, 



hj,k= 1,2,3 



O- + /3 • TT 



(12) 



(13) 



Therefore we see that the parameters a, correspond to 
a vector transformation, while f3 to an axial transforma- 
tion. This is a very important property, since this 0(4) 
symmetry can easily be generalised to 0{N) one, involv- 
ing N pion fields. 



1. Spontaneous symmetry breaking 

In the linear sigma model, the classical ground state 
is determined by the minimum of the potential term in 
the Lagrangian which describes the self-interaction of the 
scalars 
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(14) 



If the mass term /z^ is negative there is a minimum of the 
potential which is located at 



(15) 



This defines a 3-sphere, , in the 4-dimensional space 
formed by the scalar fields. Each point on is invariant 
under 0(3) rotations. For example, the point (w, 0, 0, 0) is 
invariant under the rotations of the last three components 
of the vector. This means that after a point on is 
chosen to be the classical ground state, the symmetry is 
broken spontaneously from 0(4) to 0(3) We should point 
out that the vector SU (2) isospin symmetry is isomorphic 
to 0(3) . 

Quantizing the theory we will need to expand the fields 
around the classical values, however we may choose a 
particular ground state such as 



< a >= V 



and 



< TT >= 



(16) 
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The constant quantity v is usually called the vacuum ex- 
pectation value (VEV). Then we see that we can expand 
the sigma field around the minimum as cr = w + cr', so 
the Lagrangian becomes 



C = 



1 



1 



Xv , 



(17) 



Aw , / ,2 , 2\ /2 , 2 

+ipi"/'^dfj_'>p + gvipip + igipTj^^ ■ tt + gipipcr' . 

We observe that the fermions are massive while the 
pions appear massless. This is a consequence of the 
Goldstone theorem. According this theorem spontaneous 
symmetry breaking of a continuous symmetry will result 
in massless particles, or zero energy excitations. 

After SSB, the original multiplet (cr, tt), splits into 
massless pions and a massive sigma. Also the fermions 
acquire mass which is proportional to VEV. Therefore, 
although the interaction is SU{2)l x SU{2)ii symmetric, 
the particle spectrum is only isospin SU{2) symmetric. 
This is the typical consequence of SSB. In some sense, the 
original symmetry is realized by combining the SU (2) 
multiplet with the massless Goldstone bosons to form 
the multiplets of SU{2)l x SU{2)r. The sigma field is 
massive with mass = —2fi^ — Xv^ /3 and there are 
trilinear a-ir couplings proportional to v. Another im- 
portant point is that after the SSB, the axial current will 
have a term linear in tt field, 



= ivdf^TTi H 

which is responsible for the matrix element, 

< 0\Af\TTj (p) ip^W . 



(18) 
(19) 



Using this matrix element in vr decay, we can associate the 
VEV V with the pion decay constant /tt. This coupling 
between axial current A'^ and tt^ causes the appearance 
of a massless pole. 

In the scalar self interaction, quartic, cubic, and 
quadratic terms have only 2 parameters, the coupling 
constant A and the mass term fi. This means that these 
three terms are not independent, and there is a relation 
among them. This is an example of the low energy the- 
orem for a theory with spontaneous symmetry breaking. 
We discuss briefly the low energy theorem in Section ll B 31 



2. Explicit breaking of chiral symmetry 

The SU{2)l x SU{2)fj symmetry of the linear sigma 
model is explicitly broken if the potential V{a, tt) is made 
slightly asymmetric, e.g. by the addition of the term 



A = 



(20) 



to the basic Lagrangian in Eq. To first order in the 
quantity e, this shifts the minimum of the potential to 



X 



1/2 



2^2 



(21) 




(b) 




(c) 



FIG. 1: Tree-level graphs contributing to pion scattering 
at zero temperature. Continuous lines represent pions and 
dashed lines correspond to sigma. Time runs from left to 
right. 



As a result the pions acquire mass given by 



(22) 



3. Low energy theorem 

This theorem is one of the most distinctive relations 
among amplitudes involving Goldstone bosons at low en- 
ergies. These relations are a consequence of the fact that 
Goldstone bosons are massless. Since Goldstone bosons 
do carry energies, this is possible only to the limit that 
Goldstone bosons have zero energies. Alternatively we 
can think of zero-energy Goldstone bosons as infinitesi- 
mal chiral rotation of vacuum. In a chirally symmetric 
theory, it has no effect. 

Consider the following process involving pion elastic 
scattering with or without sigma exchange. The tree- 
level contributions come from the diagrams in Fig. ^ 
The amplitudes for these diagrams are given by. 



iMa = 

iMb = 



-iX , 

-2iXv 
6 

-2iAi; 



(23) 



t — mi 



-2iXv\' 



U — TOj. 

where s, t, and u are the usual Mandelstam variables, 
s=(pi+P2)', i=(pi-P3)', U= (pi-p4)' . (24) 



Then the total amphtude is 

iM = iMa + iMb + iMc + iMa (25) 



In the hmit where pions have zero momenta, we get 
s, u — > 0, and the amphtude is 

, 3 

iM = -i\ + — . 26 

9 TOJ. 

For our choice of parameters of the hnear sigma model, 
3m^ — . Thus, the amplitude vanishes in the soft 
pion limit. 

This is equivalent to taking the limit, TOo-^ — )■ oo. Soft 
pion means that the pion momentum is much smaller 
than the sigma mass ma- These are simple examples of 
the low energy theorem which states that physical am- 
plitudes vanish at the limit where the mass of Goldstone 
bosons goes to zero. 



6 



o 8 e 




(b) (c) (d) 




(e) (f) (g) 



FIG. 2: Examples of the various types of IPI diagrams which 
contribute to the effective potential for a A(/!>^ theory. In par- 
ticular, diagrams of the type (f), are called "daisies" and those 
of type (g), super-daisies, or foam diagrams respectively. 



C. The thermal effective potential 

1. The conventional effective potential 

The finite-temperature effective potential V{4>,T), is 
defined through an effective action F((/)), which is the gen- 
erating functional of the one particle irreducible graphs (a 
graph is called one-particle irreducible (IPI) if it cannot 
become disconnected by opening only one line, otherwise 
it is one-particle reducible), having the meaning of the 
free energy of the system. Diagrammatically, the loop ex- 
pansion of the effective potential could appear as given 
in Fig. [5] Calculations using the loop expansion are very 
difficult beyond two loops. One way out of this problem 
is to perform selective summations of higher-loop graphs. 

One way to systematise such summations is the large- 
TV method, in which one considers an A'^-component field, 
and uses the fact that in the large N limit, some multi- 
loop graphs give greater contributions than others. For 
example in the A$'*, iV-component scalar theory, the 
leading multi-loop contributions come from daisy and 
superdaisy graphs in Fig. |21 




(e) (f) (g) 

FIG. 3: a-e: Examples of two-particle irreducible graphs 
which contribute to the effective potential in the CJT method, 
up to three loops. The lines represent "dressed" propagators 
G(x,y). f-g: These graphs are two-particle reducible. They 
do not contribute in the CJT method, however they do con- 
tribute to ordinary effective action In this case the lines 
represent the tree-level propagator D{(f);x,y). 



2. The Cornwall- Jackiw-Tomboulis method 

Another way to perform systematic selective summa- 
tions is to use the method of the effective action for com- 
posite operators ^2§\. In this case, the effective action is 
the generating functional of the two-particle irreducible 
(2PI) vacuum graphs (a graph is called "two-particle irre- 
ducible" if it does not become disconnected upon opening 
two lines). Now, the effective action r((/), G), depends not 
only on 4'{x), but on G{x,y), as well. These two quanti- 
ties are to be realized as the possible expectation values 



of a quantum field ^(x) and as the time ordered product 
of the field operator T^{x)^{y) respectively. This for- 
malism was initially used for a study of the 0{N) model 
at zero temperature but it has been extended to 
finite-temperature by Amelino-Camelia and Pi, and was 
used for investigations of the effective potential of the 
A(/)^ theory [2^, and gauge theories [25l| . 

There is an advantage in using the CJT method to cal- 
culate the effective potential in certain approximations 
as is, for example, the Hartree approximation of the X<f>^ 
theory. According to reference |24| , if we use an ansatz 
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for a "dressed propagator" , we need to evaluate only one 
graph, that of the "double bubble" in Fig.[3Ji., instead of 
summing the infinite class of "daisy" and "super-daisy" 
graphs, given in Figs|2f,g, using the usual tree level prop- 
agators. 

We demonstrate the advantage of this method in Sec- 
tion II El where we calculate the finite-temperature ef- 
fective potential for one scalar field, with quartic self- 
interaction. The calculation of the effective potential by 
using the CJT formalism is reviewed, in detail, in ref- 
erences [23, H^, but in order to illustrate the method 
for the calculation of the effective potential for the linear 
sigma model at finite-temperature, we will reproduce the 
basic steps here. For this presentation we follow closely 
the derivation as in [2^. 

In order to define the effective action for composite 
operators, we can follow a path analogous to the one 
leading to the ordinary effective action. The essential 
difference is that the partition function depends also on 
a bilocal source K{x,y), in addition to the local source 
J{x). As an example, we consider the theory with 
Lagrangian 



£($) 



1 



2^2 



■ TO $ 
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.$4 



(27) 



According to CJT method [23 , the generating functional 
for the Green functions in the presence of sources J{x) 
and K{x, y) is given by (we set H = 1) 



Z{J,K) = 



= / P-^exp 



2 



(28) 

where W(J, K), is the generating functional for the con- 
nected Green functions, while /(<&) — J dx'^ ^{x) is the 
classical action. Adopting a similar notation as Rischke 
and Lenaghan Q , we have used the shorthands 

(29, 

^K^ = j d'^xd*y^{x)K{x,y)^{y) . (30) 

The expectation values for the one-point function, (l){x), 
in the presence of a source J{x) is given by 



6W{J, K) _ 



5J{x) 

while the connected two-point function, G(x, y) is 
5W(J, K) _ 1 



(31) 



5K{x,y) 2 



[G(x,y)-l-0(x)0(2/)] . (32) 



The effective action for composite operators r((/), G), 
is obtained through a double Legendre transformation of 
W[J,K) = \nZ{J,K) 

r(0,G)=W(J,X)-0J-i0if0-iGif , (33) 



where GK = J d'^xd'^yG{x,y) K(y,x). Then it follows 
that 



6(j}{x) 

^^('^'g) = ^Kixy) 
SG{x,y) 2'^^'''y>- 



= -J{x)- / dSK{x,y)(^{y) , (34) 



(35) 



Physical processes correspond to vanishing sources, so 
the stationarity conditions which determine the expecta- 
tion value of the field 95 (x), and the (dressed) propagator 
Q(x,y) are given by 



5T{4>, G) 
5<t>{x) 

5n4>,G) 



5G{x,y) 



=v,G=g 



=v,G=g 



= 0, 



= 



(36) 
(37) 



As it was shown by Cornwall-Jackiw-Tomboulis in 
[23| . the effective action r(0, G) is given by 



r(0,G) = /((/.)- ^Tr (In G- 
-iTr(Z?-iG-l) 



(38) 



r2(0,G) 



where in this last equation r2(0, G) is the sum of all 
two particle irreducible (2PI) diagrams in which all lines 
represent full (dressed) propagators G, while is the 
inverse of the tree-level propagator. 



D-\x,y- 



S(j){x) S(j){y) 



(39) 



<P=v 



In the case when translation invariance is not broken, 
it can be assumed that the field (f>{x) = constant — (p. 
Then an infinite volume factor O arises from space-time 
integrations and it is customary to introduce the general- 
ized effective potential V{(j), G) = r((/), G)/Vl. The exact 
expression for the effective potential is then 



y(</),G) = f/(^) + - / InG-i(fc) 



(40) 



[D-\k;<j>)G{k)-l]+V2{<j^,G) 



where is t/((/)) is the classical potential of the Lagrangian, 
D~-^{k;(t)) is the inverse tree level propagator, which for 
A(/)^ theory is given by 



D~^{k; <P) = k^ + ^A^^ 



(41) 



and V2('/'jG) — V2{4>iG) /il. The stationarity conditions 
for the effective action reduce to the following ones in- 
volving the effective potential 



jy(0,G) 



5V{<P, G) 
6G{k) 



= 



(42) 
(43) 
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We then obtain the following relation 



(44) 



which corresponds to a Schwinger-Dyson equation for 
the full (dressed) propagator Q and where the self-energy 
S(fc) is given by 



SG{k) 



(45) 



4=<p, G=G 



The thermal effective potential has the meaning of the 
free energy density and it is related to thermodynamic 
pressure through the equation 



P 



(46) 



D. Finite temperature formalisms 

Whatever the method one chooses to calculate the ef- 
fective potential, the calculation has to be performed in 
the framework of finite temperature field theory. There 
are three main paths to follow namely: thermo-field dy- 
namics, imaginary time formalism and real time formal- 
ism. In our calculations, we do not use the machinery of 
thermo-field dynamics, but we use both real and imag- 
inary time. Each of these has certain advantages and 
disadvantages. Traditionally, the imaginary time seems 
to be convenient for systems in equilibrium, since the real 
time is thought to be the appropriate for calculations in 
systems far from equilibrium. 



1. Imaginary time formalism 

The imaginary time formalism, which is also known 
as the Matsubara formalism, provides a way of evalu- 
ating the partition function perturbatively using a di- 
agrammatic method which is analogous to that which 
is used in conventional field theory at zero tempera- 
ture 0, I23, 112, H^- According to this technique, we 
work in Euclidean space-time and use the same Feyn- 
man rules as at zero-temperature, but when evaluat- 
ing momentum space integrals, we replace integration 
over the time component ki with a summation over dis- 
crete frequencies. That means that in the case of bosons 
ki = 2TTinT, n = 0, ±1, ±2, ±3, . . .. This is encoded into 
the following relationship 



d^k 
(2^ 



/(fc) 



-y 



(2^ 



/(27rmr,k) (47) 



where (3 is the inverse temperature, (3 — l/ksT, and 
as usual Boltzmann's constant is taken to be A:^ = 1- 
For the sake of simplicity, in the following we introduc 
a subscript to denote integration and summation over 



the Matsubara frequency sums. So in what follows we 
adopt the shorthand expression 

^ E / (^^(^TTznT, k) ^ ^ f{2mnT, k) (48) 

If one is to study systems at equilibrium, the imaginary 
time formalism is adequate, but for dynamical systems 
far from equilibrium, we need to analytically continue 
back to real time which was traded in favour of the tem- 
perature. An alternative way is to work directly in the 
real time formalism from the outset. We introduce the 
basic concepts about the real time formalism in the next 
section. 



2. Real time formalism 

This formalism has a long history of applications in 
condensed matter systems, where early papers appeared 
in the 1950's. It is also known as Keldysh or closed-path 
formalism 33]. It has only been quite recently that re- 
searchers realized its usefulness in applications in particle 
physics at finite-temperature. The formalism is given in 
details in [s^. , and is also reviewed in 'sg'I , where a 
formal account for both - real time and imaginary time 
- formalisms is given. There is also a lot of information 
in the recent books 0, |^ . Further information can 
also be found in recent papers |^7, 38, 39] . 

The basic point is that for a scalar field theory, the 
propagator becomes a two by two matrix and is divided 
into two parts as 

W = Wo + iVp , (49) 
where the zero temperature part appears as 



Wo = 











(50) 



and the temperature dependent part is of the form 



where 



cosh^ e = 



2 sinh^ d sinh 29 
sinh 261 2 sinh^ 9 



=./31fco| 



(51) 



(52) 



e/3|feol _ 1 ■ 

In the limit e — > 0, the pre-factor of the thermal part 
reduces to a delta function. 

In particular the (1,1) component is given by 



+ 27r(5(fc2 - m^)n{ko) 



k^ — m? + ie 

where n(fco) is the Bose-Einstein distribution function 

1 



(53) 



n(fco) 



(54) 



exp(/3|fcol)-l ■ 
In this last equation /3 = 1 /T is the inverse temperature, 
while fco = Vk^ -f Tin? . 



E. The X^* theory at finite temperature 

As an illustration, we use the CJT formalism to cal- 
culate the effective potential for a theory. The La- 
grangian is given by 



Q 



1 



C = -{d^<f)id^'<f) ^ -m'^' - -X^\ (55) 

and in order to realize the spontaneous breaking of sym- 
metry, is considered as a negative parameter. By 
shifting the field as $ <f> + 0, the "classical potential" 
takes the form 



(56) 



and the interaction Lagrangian which describes the ver- 
tices of the shifted, theory is given by 
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(57) 



The tree-level propagator which corresponds to the 
above Lagrangian density is 



(58) 



According to CJT formalism 24J, the finite tempera- 
ture effective potential is given by 

F(0,G) = (7(0) + i / lnG-i(0;fc) 



(59) 



/ [I?-i(0;fc)G(0;fc)-l] 



+ ^2(0, G) 



where U{(f>) is the "classical potential" given by Eq. (|^ . 
and V2((/>, G) represents the infinite sum of the two- 
particle irreducible vacuum graphs. We are going to eval- 
uate the effective potential in the Hartree approximation 
which means that we only need to calculate the "dou- 
ble bubble" diagram given in Fig. |31(a). The resulting 
effective potential, is therefore. 



FIG. 4: Scliematic representation of the self-energy contri- 
bution to scalar propagator in the self-consistent Hartree ap- 
proximation. The full or "dressed" propagator appears as a 
sum of the bare propagator, and the self-energy insertion. 



In graphical terms, we could say that the self-energy is 
calculated by opening one line of each diagram contained 
in the functional of the effective potential pflL I3(tI | . In 
Hartree approximation, this is illustrated in the Fig. 0]. 

The solution Go{(j); k), of the gap equation, is inserted 
back into the expression for the effective potential, result- 
ing in a potential which is a function of (j). As is stated in 
[23], using the propagator Go(0; k) for internal lines cor- 
responds to summing all daisy and super-daisy diagrams, 
using the usual tree level propagators as in |l8l | . 

Then, we can adopt the following form for the dressed 
propagator G((/); k), 



G(0;fc) 



1 



/j2 + A-p 



(62) 



where an "effective mass" M = M{(j);k), has been in- 
troduced. The gap equation for the propagator then be- 
comes an equation for the effective mass 



2 ^ 2 Jpk^+M^ ' 



(63) 



where it is obvious, since we integrate over the loop mo- 
menta, that in this approximation, the effective mass M 
is momentum independent. 

In terms of the solution Afo('/>) of the gap Eq. (|63|l . the 
effective potential takes the form 



F(0,G) - -mV + :^A04+ / lnG-i(</>;fc) 2 24 2 jp 



+ 1 / \{k^^Tn^^\M^)G{i,-k)-\\ 



2^ ° 2^ fc2 + ^2 



/5 



A, 



1 



-t--A[ y^G(0;fc)]- 



(60) 



1 



8 ^jfi k^ + Ml 



(64) 



Minimizing the effective potential with respect to the 
"dressed propagator" G{(j);k), we obtain the gap equa- 
tion 

G"i(0; fc) = fc2 + + iA02 + / G((/.; k) . (61) 



Performing the Matsubara frequency sums as in |l8| | 
(we repeat these steps in appendix A), the logarith- 
mic integral which appears in the above expression for 
the effective potential divides into two parts. A zero- 
temperature part, Qq{M), which is divergent, and a 
nonzero-temperature part, Qp^M), which is finite, and 
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can be written as 
Q{M) 



Qo{M) + QpiM) 



(27r)3 2 
1 f d^k 



PJ (27r)= 



ln[l - exp(-/3c^k)] , (65) 



where uj^ = (k^ + M^)^/^. In this above expression, and 
in what follows, we omit the subscript on M. Similarly, 
the second integral is divided into a zero-temperature 
part Fo{AI), and a finite-temperature part Fp{M), given 
by 



F(M) 



1 



fc2 + A-P 
= Fo{M)+Fp{M) 
rf^k 1 



(27r)3 2wk 
d^k 1 



1 



(27r)3 cjk exp(/3wk) 



1 



(66) 



The second term vanishes at zero temperature, while the 
first term survives but it gives rise to divergences which 
can be carried out using appropriate renormalization pre- 
scriptions |24l l26| | . If one is interested in temperature in- 
duced effects only, as is our approximation, the divergent 
integrals can be ignored. In this case, by making a change 
to the integration variables, the finite-temperature part 
of F{AI) can be written as 



rp2 PCO 



2+2;2)l/2exp(a;2+y2)l/2_i ' 

(67) 

where we have used a shorthand notation, y ~ MjT . 

In order to simplify the lengthy expressions involved 
in what follows, we introduce the function g{z^ a) defined 
as 

" (z2 + a2)i/2 exp(z2 + a2)i/2_i ■ (68) 

Then the function /(a), will appear as the shorthand of 
the integral 



/(a) 



X dx 



1 



(69) 



(x2 + a2)i/2 exp(a;2 + a^yP - 1 
x'^ dx g{x, a) , 

and in the limit of a — 0, we recover the well known 
result /(O) = 7r2/6. 

Similarly, the finite-temperature part of the logarith- 
mic integral becomes 

n4 



Qp{M) ^—r x'^dxlnll- exp[-(.T 
27r^ Jo ^ 



y 



2U/2- 



Then, the finite-temperature effective potential, ignoring 
quantum fluctuations, can be written as 



F(0,Af) = lm2</)2 + -LA04 + Q^(Af) 

-i(M2-m2-^02)^^(M) 



(71) 



We can obtain a more compact form if we make use of 
the gap equation JB^), 



1/(^,M) = im202 + 1a<^4 + Qp{M) - ^-\[Fp{M)\ 



(70) 



(72) 



II. HARTREE APPROXIMATION 

A. Introduction 

In this section, we apply the CJT method in order to 
calculate the effective potential of the 0(4) linear sigma 
model. The a field can be used to represent the quark 
condensate, the order parameter for the chiral phase 
transition, since both exhibit the same behaviour under 
chiral transformations [Tol fisjl . The pions are very 
light particles, and can be considered approximately as 
massless Goldstone bosons. We use the model as an ef- 
fective theory for QCD, ignoring the fermion sector for 
the moment, and concentrating on the meson sector. We 
examine this model in two approximation schemes: (a) 
in the chiral limit considering the pions true Goldstone 
bosons, and therefore massless, and (b) treating them as 
massive. 



B. The linear sigma model 

As we mentioned in the introduction, the linear sigma 
model serves as a good low energy effective theory 
in order for one to gain some insight into QCD. Re- 
cently, it has attracted much attention, especially in 
studies involvin g di soriented chiral condensates (DCC's) 
El El ElEi- The model is very well suited for 
describing the physics of pions in studies of chiral sym- 
metry. Fermions may be included in the model either 
as nucleons, if one is to study nucleon interactions, or 
as quarks. The mesonic part of the model consists of 
four scalar fields, one scalar isoscalar field which is called 
the sigma a field and the usual three pion fields 7r°, tt^, 
which form a pseudo-scalar isovector. The fields form a 
four vector (ct, tt^), i = 1,2,3, which we regard as the 
chiral field and the model displays an 0(4) symmetry. 

The meson sector Lagrangian of the 0(4) sigma model 



11 



'-2 2 

—m a 
2 



^22 
— TO TT 
2 



(73) 



The last term, ecr, has been mtroduced into the above 
expression in order to generate masses for the pions. It 
is related to pion mass as £ = /tttoJ, where /tt = 93 MeV, 
is the pion decay constant. In the absence of the last term 
which breaks the chiral symmetry, the pions are massless. 

The coupling constant A of the model can be related 
to zero temperature properties of the pions, and sigma 
through the expression 



A 



3(to2 



J 7T 



(74) 



The negative mass parameter to^, is introduced in or- 
der to obtain spontaneous breaking of symmetry and its 
value is chosen to be 



-m^ = (ml -3to2)/2 > 



(75) 



To define the numerical values of A and to^, we use 
TO^ « 138 MeV, since for the sigma mass we adopt 



600 MeV. 



C. The chiral limit e = 

In order to deal with the exact chiral limit first, the 
starting point is the Lagrangian given in Eq. (|73(l . and 
we ignore the symmetry breaking term for the moment. 
The parameters A and appearing in the Lagrangian 
are simply 



A 



3to^ 



By shifting the sigma field as a ^ a 
"classical potential" is 



-2m^ ^ml>0 . (76) 
the resulting 



1 



2 i2 
TO (p 



24' 



(77) 



and the interaction Lagrangian which describes the ver- 
tices of the new theory takes the form 



A 

24^ 



24' 



A_ 

12^ 



2 2 

(T TT - 



A 

6 ' 



(78) 



In the Hartree approximation, we do not consider inter- 
actions given by the last two terms in the Lagrangian. 
We attempt to include these interactions in Section Hvl 
The tree level sigma, and pion propagators, corre- 
sponding to the above Lagrangian are 



■ TO 



(79) 
(80) 



(3) 




(9) 



(6) 



FIG. 5: Graphs of the "double-bubble" type which contribute 
to the effective potential for the 0(4) linear sigma model in 
the Hartree approximation. The numbers show the weight of 
each type of bubble, in the expression for the effective poten- 
tial. 



We evaluate the effective potential in the Hartree ap- 
proximation, which means that we only need to calculate 
the "double bubble" diagrams as in Xcf)^ theory. In the 
linear sigma model, the corresponding effective potential 
at finite temperature can be written as 



V{<p,G) 



U{<l>) + \j\nG-\cl>-k) 



- / lnG-i(0;fc) 



[V-\(l); fc)G,(0;fc) 



+ - / [p-i(</);fc)G,(0;fc) 
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4-1/2(0, G„ GO 



1] 
3] 



(81) 

where the first term U {(/)), is the classical potential and 
the last term V2((/>, Go-, G^), originates from the sum of 
"double bubble" diagrams. There are four types of dou- 
ble bubbles as we show in Fig. |S1 and these contribute 
the following terms in the potential 



V2((/), G^,G^) = 3^ 



G,(0;fc) 



+6 



24 



+I5A 
24 



/3 



G,(0;fc) / GO0;fc) 

^ 2 

G.(0;fc)' 



/3 



(82) 



Minimizing the effective potential with respect to the 
"dressed" propagators, we get the following pair of non- 
linear gap equations 



G-i(0;fc) 



T^^' + ^l G.(0;fc)-f ^ I G,{c^;k) 



5A 



G,(0;fc) 



GM: k) 



The bare propagators and are given by Eq. H79|l . 
In order to solve this system, we can use the same ansatz 
for the dressed propagators as in the one field case 



G-] =k^ + , 



(83) 
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Then by using Eqs. (|79|l . (|83|l the nonhnear system for 
the dressed propagators reduces to the foUowing system 
for the thermal effective masses 



Ml = 



\F{hU) (84a) 
^F(MO .(84b) 



In these last two equations, we have used a shorthand 
notation and F{M^^^) is given by 





(a) 



(d) 



F{M,,^) 



(85) 



As in A^"* theory, the thermal effective masses are inde- 
pendent of momentum and functions of the order param- 
eter (p, and the temperature T. 

By using these two equations, the effective potential at 
finite temperature can be written as 



y(</),M) = 



.}.(^Ml-m^-\H^)F{M,) 
-\{Ml^m^-\H^)F{AU) 



--F{M^)F{AU) 



Mi 



(86) 



Minimizing the effective potential with respect to the 
"dressed" propagators, we have found the set of non- 
linear gap equations for the effective particles' masses, 
given by Eq. H84() . In addition, by minimizing the poten- 
tial with respect to the order parameter, we obtain one 
more equation 



1 



A 



A 



= + -X<jy' + -F{M„) + -F(M,) 



(87) 



In order to study the evolution of the potential as a 
function of temperature, we perform the Matsubara fre- 
quency sums, as in the one field case. There are some 
problems concerning the renormalization of the model 
|24. lig . 1471 lig ] . At the level of our approximation, we ig- 
nore quantum fluctuations for the moment, and keep only 
the finite temperature part of the integrals. This deci- 
sion is justified by the fact that the finite terms - given 
by Eq. (|67|l and Eq. (|7(J|I - correspond to interactions 
between thermalUy excited mesons with other mesons 
present in the plasma In doing so, we neglect the 

quantum fluctuations of the meson fields and keep only 
the thermal fluctuations present in the hot plasma. The 
same choice was also adopted in [s^l ■ In a very re- 
cent investigation of the same model, and based on the 



FIG. 6: Schematic presentation of selfenergy contribution to 
pion propagator. Dashed lines represent sigma propagator, 
while continuous ones correspond to pions. 



(«) 



(&) 



Q C)^ 



id) 



FIG. 7: Schematic presentation of selfenergy contribution to 
sigma propagator. Dashed lines represent sigma propagator, 
while continuous ones correspond to pions. 



same formalism, Rischke and Lenaghan [51^ have shown 
that the model is indeed renormalisable. Renormaliza- 
tion of the model has also demonstrated by Chiku and 
Hatsuda |5^|5i|, who on the basis of optimised perturba- 
tion theory, and the real time formalism, have studied the 
spectral densities and the effective masses of the sigma 
and the pions. 

Diagrammatically, differentiation of the double bub- 
bles in Fig. O gives the selfenergy loop contribution to 
the particle propagators as we show in Figs. 16171 We can 
picture the dressed propagator with a thick line which ap- 
pears as a sum of the bare propagator (not shown in the 
figures) and the selfenergy contributions. Actually one 
can see these two figures as a graphical representation 
of the gap equations for the dressed propagators. As we 
have mentioned already, at finite temperature and in the 
real time formalism, the propagators appear to consist of 
two parts, a zero temperature part and a thermal part. 
We denote each zero tempearture loop with continuous 
line and each thermal with a cut. Our choice to keep only 
the non divergent parts of the integrals, is equivalent to 
taking into account only the thermal propagators. This 
means that for pions, we are taking into account only the 
thermal self energy loops in Figs. EId, d while for sigma 
those in Figs. [3),d. 
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Using a compact notation, the finite temperature ef- 
fective potential can be written in the form 

V{c^,M) = lmV + ^A04 

-^F^{M,)Fi3{M^) , (88) 

where in this last step, we have used the gap equations 
given by Eq. 1)84(1 . The exact expressions for Fp{M) and 
Qp{M), are given by Eq. and Eq. ifTUI) respectively. 



1. High temperature limit 

In order to calculate the effective masses as functions 
of temperature, we need to solve the system of the three 
equations given by Eq. and Eq. l|57|) . We first ob- 
serve that if = 0, which happens in the high temper- 
ature phase, the two equations become degenerate, the 
particles have the same mass, and we have to solve only 
one equation 



(89) 



As in the expression of the effective potential, we keep 
only the finite temperature part of the integral. This last 
equation can be used to define a "transition temperature" 
Tel. This temperature is defined as the temperature, 
where both particles become massless. Recall now that 
Ffj{M) is given by 



Fp{M) 



27r2 



dx 



1 



ai/2 exp(ai/2) _ i 



(90) 



where a = -|- j/^ and y = M/T . When the mass of the 
particles vanishes, this integral reduces to the well known 
result 



I{x) = 



■ dx 



~6 



(91) 



Therefore by using Eqs. (jSHl), 1101) and we find that 



= TO^ + A 



27r2 6 



(92) 



Actually this defines the transition temperature T — Td 
as 



T,i = V2 - 



6m^ 



2\ 1/2 



(93) 



But, in defining our model parameters, we have chosen 
that at zero temperature (f)^ = f'^ = —Qrn^/X, where 
/tt = 93 MeV is the pion decay constant, so we find that 
Tci = V2/^ w 131.5 MeV. 
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FIG. 8: Solution of the system of gap equations in the chiral 
limit, e = 0. The sigma and pion effective masses are given 
as functions of temperature. 



2. Low temperature limit 

In the low temperature phase, we can eliminate cf) be- 
tween Eqs. (|H3I) and ^J^. We then end up with the 
following nonlinear system 



Ml = -2m^ - XFpiAU) - XFpiA'U) 



Mi 



A 



;FpiM, 



A 



Fp{M^) 



(94) 



where, as before, we have ignored the divergent parts of 
F[Mfj/T^). This system has been solved numerically using 
a Newton-Raphson method, and the solution is presented 
in Fig. IHl . 

As shown in Fig. |S1 the temperature Td which is calcu- 
lated numerically, was found to be in excellent agreement 
with the value obtained by using the limit of the high 
temperature equations with degenerate masses. 

At this point, we can observe that there is an indica- 
tion of a first order phase transition, because combining 
Eq. l|84H|) with Eq. we find that 



Mi 



(95) 



This last equation shows, of course, that the order pa- 
rameter varies with temperature proportionally to the 
sigma mass. The temperature dependence oi (j) = (f>{T), 
is calculated by using the sigma mass as it was found by 
solving the system in Eq. (|94|) . This is shown in Fig. |51 
where it is obvious that this approximation predicts a 
first order phase transition, because (f>, which is the or- 
der parameter of the phase transition, appears to have 
two different values for the same temperature. 

This last observation coincides with the qualitative pic- 
ture given by Baym and Grinstein, in their early paper 
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FIG. 9: Evolution of the order parameter as a function of 
temperature. 



[43, where their "modified Hartree approximation" pre- 
dicts a first order phase transition, as well. However, in 
contrast to our approximation, they had included quan- 
tum fluctuations into their analysis as well. Signals of 
a first order phase transition have also been reported in 
recent analyses by Randrup Roh and Matsui [i^ . 
and Rischke and Lenaghan [5l[. However, this result 
seems to disagree with other investigations of the linear 
sigma model (including fermions) or the Nambu-Jona- 
Lasinio model, where a second-order transition has been 
reported ;52. .5.1. .54i| . 



D. Evolution of the thermal potential 

In order to get more insight into the nature of the phase 
transition and verify that the transition is of the first or- 
der, we can calculate the effective potential V{(j),T) as a 
function of the temperature and the order parameter. So, 
we first solve numerically the system of three equations 
in Eq. (|84|l and Eq. (|87|l (where of course we keep only 
the finite temperature part of the integrals), and calcu- 
late the effective masses of the particles as functions of 
the order parameter and the temperature. Finally, the 
effective potential is calculated numerically, using these 
masses. The evolution of the potential for several tem- 
peratures is given in Fig. ^| The shape of the potential 
confirms that a first-order phase transition takes place, 
since it exhibits two degenerate minima at a temperature 
Tc « 182 MeV, which is usually defined as the transition 
temperature. The second minimum of the potential at 
4> disappears at a temperature Tc2 ~ 187 MeV . 

Normally, in first order phase transitions there are also 
two other temperatures of interest, apart from the tran- 
sition temperature, but the relevant isotherms are not 
shown explicitly in this picture. These temperatures 
Tci and Tc2 are called in condensed matter terminology, 
the lower and upper spinodal points, respectively. Be- 
tween these temperatures, metastable states exist, and 
the system can exhibit supercooling or superheating. For 
Tci < T < Tc, the metastable states are centred around 




FIG. 10: A qualitative picture of the evolution of the effective 
potential V{cf>, T) as a function of the order parameter for 
several temperatures. For our choice of parameters of the 
linear sigma model, the two minima appear as degenerate at 
Tc « 182 MeV, while the upper spinodal point (not shown 
explicitly) appears at Tcj ~ 187 MeV . 



the origin, since for Tc < T < Tc2 the metastable states 
occur for 7^ 0. When the system reaches Td or Tc2, the 
curvature of the potential at the metastable minima van- 
ishes. A discussion about first-order phase transitions, 
and more details_a_bout how these transitions proceed, 
can be found in 



E. The broken symmetry case e 7^ 

When e 7^ 0, the term linear in the sigma field into the 
Lagrangian generates the pion observed masses. This 
term is G independent and so, minimization of the po- 
tential with respect to "dressed" propagators will give us 
the same set of gap equations for the effective masses as 
before. However, minimizing the potential with respect 
to (p, we get the following equation 



i-e = . (96) 
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FIG. 11: Solution of the system of gap equations in the case 
when e 7^ 0. At low temperatures the pions appear with the 
observed masses TOtt ~ 138 MeV . 



In order to proceed, we need to solve the nonlinear system 
of three equations in Eq. and Eq. (|^ . We first 

observe that at T = 0, Eq. H96|l becomes 

Ml = m^ + \\<i>^ = t (97) 
6 

where m,r is the tree level pion mass. Then for = 
we recover the relation between the pion mass at zero 
temperature and the symmetry breaking factor e: e = 
J-n'rn^, where /tt is the pion decay constant. We solved 
the system of Eq. (|84|l and Eq. (|96|l numerically, and the 
solution is presented in Fig. ^2 . 

At low temperatures, the pions appear with the ob- 
served masses, but their mass increases with tempera- 
ture since the sigma mass decreases. At high temper- 
atures (higher than ~ 300 MeV), due to interactions in 
the thermal bath, all particles appear to have the same 
effective mass. 

The presence of the symmetry-breaking term into the 
Lagrangian, modifies the evolution of the order parame- 
ter </), as well. As it is obvious in Fig.^l as the tempera- 
ture increases, the order parameter decreases, and at very 
high temperatures vanishes smoothly. But in this case, 
the change is not a phase transition any more. We rather 
encounter a smooth crossover from a low-temperature 
phase, where the particles appear with different masses, 
to a high-temperature phase, where the thermal contri- 
bution to the effective masses makes them degenerate. 



III. LARGE TV APPROXIMATION 

A. Introduction 

The large- A'' approximation of the linear sigma model 
has been studied recently by Amelino-Camelia , and 
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FIG. 12: Evolution of the order parameter </!> as a function of 
temperature. 



our expressions are very similar to the ones obtained 
there, since the same method is used in both cases. How- 
ever, in our approach, we do not consider the renormal- 
ization of the model because in our approximation we 
take into account only finite temperature effects. Our 
analysis, is in a sense, complementary to that in |46j| . 
since we solve the system of gap equations and consider 
the effects of the symmetry breaking term (the last term 
in the Lagrangian given by Eq. (|98|l 'l. which is omitted in 
reference |46j . The renormalization of the model is inves- 
tigated in a recent publication by Rischke and Lenaghan 
|51| . Another recent treatment of the 0{N) model ap- 
pears in [t^, where they use the CJT method, as well, 
in order to calculate the thermal effective potential. The 
0{N) version of the linear sigma model is very popular 
in condensed matter studies and it is a very well-studied 
model. Other studies of this model in particle physics 
context are those in [63l l6^ . 

The generalised version of the meson sector of the lin- 
ear sigma model is called the 0{N) or vector model, and 
is based on a set of N real scalar fields. The 0{N) model 
Lagrangian can be written as 



1 



1 



-TO f> 



1 



and in the absence of the last term, it remains invariant 
under 0{N) symmetry transformations for any N x N 
orthogonal matrix. Our choice of taking negative, 
results in a A'^-dimensional "mexican hat" potential. 

In order for our notation to be consistent with applica- 
tions to pion phenomenology, we can identify <I>i with the 
a field and the remaining iV — 1 components as the pion 
fields, that is $ = (c, tti, . . . ,7rjv_i). The last term, ea 
in the above expression has been introduced in order to 
generate masses for the pions. When = 4, this model 
is exactly the linear sigma model of the previous section. 

We proceed by examining the 0{N) within two ap- 
proximation schemes, as we did in the Hartree approach. 



16 



in the previous section. 



temperature as 



B. The chiral limit e = 

In order to study the 0{N) model, we proceed in abso- 
lutely analogous steps as in the Hartree case in SectionlTll 
So, by shifting the sigma field as a ^ a + (j), the tree level 
propagators are 

2A 



D-\<j,;k) = e+m^ + —cj,' 
2A 



(99) 



Then the effective potential at finite temperature will 
appear as 



M) 



1 



2 j,2 

m (p 



1 

m 



lnG;i(0;fc) 



N~ 1 



lnG-i(0;fc) 



i J [2?-i(0;ft)G.(</>;fc) 



N -1 



[p-i(0;fc)G,(0;fc)-(7V-l)] 

(100) 



+ V2('/>, Go-, Gtt) , 



where the last term originates from the double bubble 
diagrams and its contribution is 



F2(0,G,,G,) 



A 

67V 



Go(0;fc) 
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+ 



XjN^ - 1) 

6iV 
A(iV- 1) 
6^V 



1 2 



G,(0;fc) 



/3 



Go(<^;fc) / G,(0;fc) 

J/3 



The weight factors appearing in the above expression can 
be understood in a similar way as in the 0(4) case, with 
the only difference being the A'' — 1 pion fields. Of course, 
it is easy to see that we recover the previous case by 
simply substituting N = A. 

As in the case of A(/)^ and the 0(4) model, we minimise 
the effective potential with respect to the dressed propa- 
gators, and we get a set of gap equations. By using the 
same form for the dressed propagators as before, we end 
up with the following set of nonlinear gap equations for 
the thermal effective particle masses 



Mi 



Mi 



2A 

37V 



2X{N + 1) 
37V 



FpiM^) , 



where we only keep the finite temperature part of the 
integrals, and we have defined the particle masses at zero 



iV^ 
2X_ 

3iv' 



(102) 



As it is easy to observe, for = 4 we obtain identical 
expressions for the system of gap equations, as in the case 
of the 0(4) model. 

In the large-A^ approximation, which means that we 
ignore terms of 0(1/A^), the system of the last two equa- 
tions reduces to 



M2 

Mi 



^0^ + ^F,(M.) 



2A 

37V' 



2A 



FfiiM^) 



(103) 



We have retained the terms quadratic in (j) since (f) de- 
pends on 7V as = — 37Vm^/2A, and so these terms are 
of 0(1). In order to solve this system, and be in "some 
contact" with phenomenology in the chiral limit, we can 
set iV = 4. Then, the pious are massless and the sigma 
has a mass M^ — ~2m^ at zero temperature. Now our 
system is written as 



7V/2 
Mi 



1 



2A 

T 

2A 



Fp{M,) 
Fp{M^) 



(104a) 
(104b) 



The effective potential will appear in the form 



1/(0, M) = im202 + J_Xcj)i + i / in(fc2 + Ml) 



2 ' 6N 
(N-l) 



[ In(fc2 + Af2 



^^iM'^-m'-^^^)F,{M^) 
2N 

X{N-1) 



3N 



67V 

Fp{M„)Fp{M^) 



(105) 



In order to solve the system for the thermal effective 
masses in Eq. (|104|l . we proceed as in the Hartree ap- 
proximation. At very high temperatures, the potential 
has only one minimum, that at </) = 0, and in this case, 
the two equations become degenerate 



Mi = Mi = M^ 



2A 

™' + ^Ffi{M) . (106) 



This last equation actually defines the critical tempera- 
ture. FplM) is given by the same expression, as in the 
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0(4) case. The mass of the particles vanishes at the criti- 
cal temperature, so we can use the result given in Eq. H9U|) 
to find that the critical temperature is at 



Tc = V3 



6m^ 



1/2 



V3U ~ 161 MeV . (107) 



Before proceeding to examine the low-temperature 
phase, we should make an observation which actually ex- 
poses the significant difference between the Hartree ap- 
proximation in the A'^ = 4 case, and the large- approxi- 
mation. Minimizing the potential with respect to gives 



dV{(f), M) 



2A , , 2A(iV-l) , , 



= , 



which, in the large-TV approximation becomes 



3 



= . (108) 



Combining this last equation with Eq. 
observe that 



104b|) above, we 



dV{(l),M) 



4>Mi = . 



(109) 



Therefore, the large-A^ approximation implies that the 
pions should be massless, in the low-temperature phase 
in accordance with the Goldstone theorem. 

This observation is reflected in the solution of the sys- 
tem of the gap equations, as is shown in Fig. ^| The 
pions at low temperatures appear as massless, but at 
high temperatures the thermal contribution to the effec- 
tive masses make them degenerate with the sigma. The 
order parameter vanishes continuously in this case as is 
shown in Fig. 1141 and this corresponds to a second-order 
phase transition. 
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FIG. 13: Solution of the system of gap equations in the large- 
approximation in the chiral limit. At low temperatures, the 
pions appear as massless. 
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FIG. 14: Evolution of the order parameter with temperature 
in the large- A'' approximation in the chiral limit £ = 



C. The broken symmetry case e ^ 

As already mentioned for the 0(4) case, the symme- 
try breaking term ea has been introduced into the La- 
grangian in order to generate the observed masses of the 
pions. The same can be done for the 0{N) model - the 
only difference being the A^ — 1 pion fields. Inserting this 
term into the expression for the effective potential and, 
differentiating with respect to we obtain, as in the 0(4) 
case, one more equation. In the large-A^ approximation 
this is written as 



e = 



(110) 



We have solved this last system of three equations 
given by Eq. H104|) and Eq. (|110|) numerically, and the 



solution is given in Fig.^|. As in the A^ = 4 case, there 
is no longer any phase transition. We encounter again 
the crossover phenomenon between the low- and high- 
temperature phases, the difference now being that the 
change of the order parameter fFig. ll6|) in the transition 
region is much smoother than the "sharper" behaviour 
seen in the A^ = 4 case, in Fig. ^] . 

At this point, we would like to comment about the re- 
sults presented in the previous two sections. Part of this 
work has also been presented elsewhere j6^ |6^ . Firstly, 
the CJT formalism of composite operators proved to be 
very handy because we actually needed to calculate only 
one type of diagram. In both cases, we solved the system 
of gap equations numerically, and found the evolution 
with temperature of the effective thermal masses. In the 
Hartree approximation, we find a first-order phase tran- 
sition but, in contrast, the large- A'^ approximation pre- 
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FIG. 15: Solution of the system of gap equations in the large- 
approximation in the case of broken chiral symmetry e ^ 
0. At low temperatures, the pions appear with the observed 
masses. 
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FIG. 16: Evolution of the order parameter as a function of 
temperature. 



diets a second-order phase transition. This last observa- 
tion seems to be in agreement with different approaches 
to the chiral phase transition, based on the argument 
that the linear sigma model belongs in the same univer- 
sality class as other models which are known to exhibit 
second-order phase transitions Q • The same conclusion 
appears to be in the work of Bochkarev and Kapusta [63j | , 
where the linear sigma model examined in the large-A'^ 
approximation, and they report second-order phase tran- 
sition, as well. 

However, in the large-A'^ approximation, the sigma 
contribution is ignored and this of course introduces er- 
rors when we calculate the critical temperature. In the 
case iV = 4, which is closer to phenomenology, we could 
have probably obtained a better approximation if we had 
considered the effects of interactions given by the last two 




FIG. 17: Solution of the system of gap equations in the large- 
approximation in the case of broken chiral symmetry for 
various values of e. We observe in this graph that as e — » 
the solution approaches the one given in Fig. ^] . 
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FIG. 18: The behaviour of the order parameter <j> as a. func- 
tion of temperature in the case when chiral symmetry is bro- 
ken (e 7^ 0) in Hartree (continuous line) and large-N (dashed 
line) . 



terms in the Lagrangian given by Eq. 1)78(1 . We attempt 
to study the effects of these terms in Section Hvl 

When we include the symmetry breaking term ea 
which generates the pion observed masses, both in 
Hartree and large-iV approximations, we found that 
there is no longer any phase transition. Instead, we ob- 
serve a crossover phenomenon where the change of the 
order parameter in the Hartree case occurs more rapidly 
in contrast to the smoother behaviour exhibited in the 
large-A'^ approximation. The difference in the behaviour 
of the order parameter in these two cases is shown in 
Fig.UHl. 



19 



This observation confirms results reported elsewhere as 
for example in the report by Smilga or in the recent 
papers by Chiku and Hatsuda pni l5q . In the latter 
analysis they also report indication of a first-order phase 
transition in the chiral limit. In a recent investigation 
of the linear sigma model by Rischke and Lenaghan , 
the results obtained are similar to ours, but they examine 
renormalization of the model as well. 



IV. BEYOND THE HARTREE 
APPROXIMATION 

A. The sunset diagrams 

In order to understand the nature of the chiral phase 
transition, in Section^ we have calculated the thermal 
effective potential in the Hartree approximation, which 
means that we did not take into account all the inter- 
actions stemming from the Lagrangian of the model. 
We found that this approximation predicts a first-order 
phase transition. However, this partial resummation of 
the Hartree approximation seems to conflict with other 
approaches, as for example, the one of Rajagopal and 
Wilczek '6','4l,f57|, where the same model is investigated 
on the basis of universality arguments and a second-order 
phase transition is found. In order to reach deeper insight 
into the nature of the chiral phase transition, we need 
take into account all the interactions that the Lagrangian 
of the model describes, at least at the two-loop level cal- 
culation of the thermal effective potential. The study of 
the effective potential beyond leading order for the elec- 
troweak phase transition has been studied by Arnold and 
Espinosa |60j while the QCD case is discussed in [6lll6^ . 

As it is shown in Section by shifting the sigma field 
as a — > a + (j), the interaction Lagrangian is given by 
Eq. (|7S|l . This interaction part describes proceses of the 

form TTTT — !■ TTTT, I7(T —^ and (TtT — > TTTT Or CTTT fJTT. 

We have considered the interactions of this form involving 
thermal pions and sigmas, when we calculated the bubble 
diagrams in the Section^] However, the last two terms in 
Eq. (|78|l give rise to processes of the form aa ^ a ^ aa 
and also to tttt — s- ct — s- tttt as well as interactions between 
pions with a sigma exchange. 

In calculating the potential at the two-loop level, one 
now is dealing with the so-called "sunset diagram". 
There are two sunset diagrams, and we show them in 
Fig. El The thermal effective potential now contains 
terms of the Hartree approximation (the double bubble 
diagrams) as in Section m plus the new terms of the 
sunset diagrams and it has the form 
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FIG. 19: The two sunset-type diagrams which contribute to 
the effective potential of the linear sigma model at two loop 
level. 
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FIG. 20: Schematic representation of the selfenergy contri- 
bution to pion propagator. Continuous lines correspond to 
pion propagator, and dashed lines to sigma. The thick line 
represents the full propagator. 

B. The gap equations 

As we have already pointed out, functional minimiza- 
tion of the effective potential with respect to the "dressed 
propagator" will give us a set of equations for the effective 
particle masses. Graphically, variation of the effective 
potential with respect to the dressed propagator, corre- 
sponds to opening one propagator line in all diagrams of 
the effective potential. Taking the functional derivatives 
of the effective potential with respect the pion full prop- 
agator Gtt, we find the gap equation for the pions. We 
can represent this equation diagrammatically, as shown 
in Fig. 1201 . Repeating this procedure for the sigma full 
propagator Ga, we find the sigma gap equation. The di- 
agrammatic representation for this gap equation is given 
in Fig. El . 

As discussed in the Section in the real time formal- 
ism, the thermal boson propagator becomes a two by two 
matrix, but at the level at which we arc working we need 
only the (1^1) component (we ommit the index 11 for 
simplicity) [Sa . l59j| . This is precisely the real time prop- 
agator, as given by Dolan and Jackiw This propa- 
gator consists of a sum of two parts: a zero temperature 
part which corresponds to 



1/2(0,T) =I/Hartrcc + 14unsct . (HI) I?o(fc) = T^ ^ , , ^ , + 7r,5(fc^ - m^) , 

(112) 
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FIG. 21: Schematic representation of the self-energy contri- 
bution to sigma propagator. Continuous lines correspond to 
pion propagator, and dashed lines to sigma. Thick line rep- 
resents the full propagator. 



where P means the principal value, and a thermal part 
given by 



VT{k,T) = 27r(5(fc2 - m'^)n{ko) 



(113) 



where n^ko) is the Bose-Einstein distribution function 



n(fco) 



exp(/3 |fco|)-l 



(114) 



where the factror /3 — l/T defines the inverse tempera- 
ture while ko — \/k^ -I- is the energy. 

The first part of the propagator is the usual virtual par- 
ticle exchange, present for all four-momenta k. The sec- 
ond term, the thermal part, of the propagator describes 
real (on shell) particles existing in the hot plasma. These 
particle particles are present only when fc^ = as the 
delta fuction constrains what follows, we denote 

the zero temperature part with continuous lines, and the 
thermal part with lines incorporating a cut " | " . 

In this part of the work, we cannot follow blindly the 
same recipe as in the Hartree case using a dressed propa- 
gator. The reason is that now the self-energy is momen- 
tum dependent as it is obvious from the graph in Fig. l2Ub . 
and the ones in Fig. I21b.d . We can avoid this difficulty, 
if we adopt a "Hartree like" recipe and accept a form of 
the dressed propagator as 



1 



M2, 
tt/o 



(115) 



where we will consider this effective mass as obtained 
from the self-energy for zero external momenta. 

If we try to distinguish the thermal from the quan- 
tum fluctuations, the pion full propagator could be rep- 
resented by the following set of graphs as is shown in 
Fig. 1221 . On the other hand the full sigma propagator 
will be represented as is shown in the Fig. 1231 . 
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FIG. 22; Schematic representation of the self-energy contri- 
bution to pion propagator. Continuous lines represent pions, 
dashed lines are sigmas, and lines with vertical dashes corre- 
spond to the thermal propagators. 
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FIG. 23: Schematic representation of the self-energy con- 
tribution to sigma propagator. Continuous lines represent 
pions, dashed lines are sigmas, and lines with vertical dashes 
correspond to thermal propagators. 
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(a) 






FIG. 24: Schematic representation of the sunset graphs. Con- 
tinuous lines represent pions, dashed lines are sigmas, and 
lines with vertical dashes correspond to thermal propagators. 



C. The pion gap equation 






(b) 





FIG. 25: (a) The sunset diagram with two thermal pions. 
A hne with a dash corresponds to the thermal part of the 
full propagator, (b)-(d) The equivalent scattering diagrams 
involving thermal pions. 



In our attempt to extend the Hartree approximation, 
our guide is the result of the low-energy theorem which 
was presented in the first chapter. We have shown that 
the invariant amplitude of the sum of all four diagrams 
in Fig. n vanishes. This acts as motivation to consider 
the equivalent in the case where we deal with thermal 
particles. As a first step, we concetrate at the low- 
temperature region and mainly on the effects of the sun- 
set diagrams on the effective masses of the pions. 

We have shown the sunset diagrams which contribute 
to the effective potential in Fig. E| We can redraw these 
diagrams indicating the thermal propagators with a cut 
as we did in the graphic representation of the gap equa- 
tions. These graphs should have the same topology as in 
Fig- El but with one, two or three cut lines representing 
thermal propagators. There is no graph with 3 thermal 
sigmas, because of energy conservation at the vertex. We 
show the complete set of sunset diagrams in Fig. |2] ■ 

As a first approximation in the resummation of the 
subclass of the sunset diagrams, we only consider the 
sunset diagram with two thermal pions. Our motivation 
stems from the fact that at zero temperature and in the 
exact chiral limit, the pions being true Goldstone bosons 
have vanishing scattering amplitude as was presented in 
the first chapter. As is shown in Fig. [23 the sunset di- 
agram with two thermal pions corresponds to the case 
where two real (thermal) pions collide to form a sigma, 
which afterwards decays into two pions. It also corre- 
sponds to pion clastic scattering with sigma exchange. 

The sunset diagram with two thermal pions is given 
in Fig. 125b . and contributes the following term in the 



effective potential 



Vsunset Ga: Gyr) 



(116) 



where Q is an shorthand form of the integral 
^7r/o-(0;fc)= / G,,/a{(f>;k) 



(117) 



Propagators indicated as G^ correspond to the thermal 
part of the real time propagator given by Eq. (|113|l . 

The full expession corresponding to the diagram in 
Fig. Oil is 



,(M,,M,) = 



5(^2 _ 



(27r)4 



27r n(ko) 2tt n(po) 



Ml 



(118) 



Following Arnold and Espinosa 60] , we adopt the conve- 
nient notation = (fco,k), for the four-momenta and, 
reserve the symbol k = |k|, for the magnitude of the three 
momentum. In this section we have adopted a different 
formalism than in Sections Hll and ITTll so the superscripts 
denote to which particle's self-energy we are reffered to, 
while the subscripts denote which particles are running 
in the loop. A subscript 6 denotes a thermal propagator. 

Taking the derivatives of the potential with respect to 
the full propagator G^, we get the diagrammatic equa- 
tion for the pion prapagator given in Fig. 1221 . As a first 
approximation we make a selective summation by includ- 
ing the two graphs as in Hartree CFigs.l^. d), plus the 
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FIG. 26: The pion self-energy diagram with one thermal 
pion. 



graph with one thermal pion fFig.l22F). The topology of 
this last graph is given in Fig. [211 • 

If we use the real time form of the propagator as in 
Eqs. H112|l and (|113|l . the graph in Fig. contributes 
with the following expression 



d'^K 2'Kn{E^)5{K^ - Ml) 



(27r)4 [P + X)2 - Af2 + it 



(119) 



Then using Eq. 11121) . we can split the above integral 
into real and imaginary parts. The imaginary part is re- 
lated to dissipation phenomena that occur as the parti- 
cles propagate within the thermal plasma. We comment 
on these matters in Section Ivl where we compare our ap- 
proach to the chiral phase transition to the work of other 
investigators. 

At the moment, we are interested only in the real part 
of the above integral, which has the form 



P(M,,M,) 



(27r)4 [P + KY - M2 



(120) 



Evaluation of this integral is given in Appendix^ Then 
the resulting gap equation for the pions is 



Mi 



1 \ 'i A 

m2 -I- + ^F{M,) + ^-F{M^) 
6 6 6 



P{M^,M,) 
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Taking the derivative of potential with respect to the 
order parameter results in the following equation 



= + \xcf>'' + ^F{M„) + ^F{M^) 
6 2 2 



A2 



R{M^,M,) 



(122) 



where the last term comes from the sunset diagram given 
in Fig. . 

Combining the above equations, we find that the ther- 
mal pion mass is given by 



A 



A 



M'^ = -^F{M,) + -F{M,) 



3 

A2 02 



A2 



g P(M,, M,) + —R{M^, M,) . (123) 



As a first step we are interested in the low-temperature 
phase, and since sigma is very heavy, we do not expect 
significant changes in the sigma mass, so we can use that 
M^ — ^XcfP. Also the contribution from F{Mcr) can be 
neglected as a first approximation, since it is exponen- 
tially supressed. Therefore in order to find out how the 
pion efective mass evolves with temperature, we need to 
calculate the terms R and P. 

Evaluating the integral which corresponds to the pion 
self-energy graph in Fig. [211 we find 



27r^ Jo K 
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Details of this calculation can be found in Appendix ^ 
On the other hand, evaluation of the sunset integral in 
Fig.[23i, results m 
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where 



X{k,p,M^,M,) 



X j pdpg{p,M^) 

X ^^^SpP:^^ (125) 
Y{k,p,M^,M^) ^ ' 



Ml{e + p2) + MliMl - Mil A) 
+ (2Af2 - Ml)kp , 



Y{k,p, M^, M„) = Ml{e + p") + Ml{Ml - Mil A) 
-{2Ml-Ml)kp 

and we have used the abbreviatiated form for g{p, M) 
introduced by Eq. (|68|l . The exact calculation is given in 
Appendix IF II . 

It is interesting to see how these integrals behave in the 
limit of vanishing pion mass. Using a compact notation 
(introduced in Sectional, F{M) can be written as 



T 



F{M,/.) = ^f{M,/jT) . 



(126) 



Rescaling the integral as a; = j3k, and in the limit of van- 
ishing pion masses, we obtain the following expressions 
for the real part of the above integral, corresponding to 
self-energy graph 
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FIG. 27: The low-temperature pion mass M^, as a func- 
tion of temperature. The continuous hne corresponds to the 
Hartree approximation. 

On the other hand, the integral corresponding to sunset 
graph is given by 

R{M^,M„) = x^dxgix,MjT) 

X / y^dyg{y,AU/T) . (128) 
Jo 

Recall now that the sigma mass is given by 

AC = iA^2 (129) 

Inserting these into the pion gap equation, we can observe 
that F(M7r) and -P(M^, Ma) cancel exactly, so we end up 
with an expession where the thermal contribution to the 
pion mass squared is ©(T^). The same result under a 
different approach has already been reported by Itoyama 
and Mueller ^ . 

We illustrate the situation with in Fig. 1271 where we 
plot the pion mass as we have calculated above, with 
the effective pion mass as was calculated in Hartree ap- 
proximation. As we can observe in Fig. 1271 in the Hartree 
calculation the pion mass is proportional to temperature. 
In contrast, inclusion of the self-energy graph results in 
the mass being ~ and obviously is small at low tem- 
peratures. 

D. The sigma gap equation 

In our attempts to consider all the self-energy graphs 
with one thermal particle in the loop, we faced serious 
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FIG. 28: The sunset diagram with two thermal sigmas. A 
line with a dash corresponds to the thermal part of the full 
propagator. 
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FIG. 29: The sigma self-energy diagram, with one thermal 
sigma. 



numerical difficulties mainly from the pion self-energy 
graph. Being unable to solve the complete system of 
equations, we have tried to understand the effects of each 
one graph individually. 

If we add to the potential the sunset with two thermal 
sigmas, which is given in the Fig. 1281 differentiating, we 
find a gap equation for sigma which contains the contri- 
bution of the sigma self-energy graph given in Fig. 1^ 

We proceed, using the same system of equations as in 
Hartree, with the only difference being that the sigma 
gap equation will contain the contribution of this graph. 
We have also introduced the symmetry breaking term, 
so the pious are massive. The result, given in Fig. 1201 is 
somehow the one expected. We know that the sigma is 
heavy so, at low temperatures, there is no significant dif- 
ference with the Hartree calculation. On the other hand, 
at high temperatures the masses will approach ideal gas 
behaviour. Only in the intermediate region of temper- 
atures we can observe a small deviation of the Hartree 
result. 

There is also no problem in calculating the sigma self- 
energy graph with one thermal pion given in Fig. 1311 
However, the addition of the relevant term into the 
sigma gap equation deviates the sigma mass from ma- = 
600 MeV as it is clear in Fig. ■ We suspect that this 
is due to the fact that we do not include the vacuum 
graphs. 
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FIG. 30: Solution of the system of gap equations for the sigma 
and pion effective masses. Dashed hnes correspond to the 
Hartree solution as in Section ini The pion gap equation is 
as in the Hartree case, while in the equation for the sigma we 
have included the sigma self-energy graph with one thermal 
sigma. 
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FIG. 31: The sigma self-energy diagram with one thermal 
pion. 



V. OTHER APPROACHES 
A. Propagation of pions in hot plasmas 

Self-consistent approximations in many body systems 
have been tried a long time ago. To our knowledge, these 
attempts date back to the 1960's, with early papers such 
as the ones by Luttinger and Ward |23|, or the ones by 
Baym 28]. Recently these ideas have been employed by 
many people in various contexts as in the recent works 
by Knoll, Van Hees and others 29, 30]. 

Our approach to the chiral phase transition is not, of 
course, the only one. There have been a lot of studies on 
the subject since it is related to important phenomena as 
the formation of quark-gluon plasma, the dilepton em- 
mision etc. People have used various ways to approach 
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FIG. 32: Solution of the system of gap equations for the sigma 
and pion effective masses. Dashed lines correspond to Hartree 
solution as in Section im The pion gap equation is as in the 
Hartree case, while in the equation for the sigma we have 
included the sigma self-energy graph with one thermal pion. 

these problems, and we are going to outline their success 
here. 

A basic ingredient for the understanding of several 
physical processes that take place in hadronic plasmas, 
is the propagation properties of pions. The properties 
of pions in hot hadronic matter are encoded in the pion 
propagator. The real part of the pion self energy is re- 
lated to dispersion and group velocity, while the imag- 
inary part encodes the information about the pion ab- 
sorption [H^) IzO; In general, the expression which 
gives the pion mass shift is of the form 

+ + , (130) 

where E is the pion self energy, and it depends on the 
physical conditions of the medium, in which the pion 
propagates. 

At the one loop level as it is the Hartree approximation 
there is no imaginary part in the self energy. The effective 
mass of the particles acquires just thermal contribution. 
However, when we include the sunset diagrams into the 
potential the mass gap equations contain terms which are 
momentum dependent. 

A recent work, which uses the linear sigma model in 
order to study dissipation phenomena in hot matter, is 
the one by Ayala, Sahu and Napsuciale [llllllzl ■ An- 
other work about dissipation properies of the pions, is the 
one by Rischke [t^. The pion self-energy is calculated 
at the one-loop level. In order to perform the study, 
the momentum running in the loop is divided into soft 
and hard. The model is used to study disoriented chiral 
condensates. 
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B. Other approaches using the sigma model 

The linear sigma model is a rich theory in its own right, 
however it has been used as an effective theory to QCD 
transition since direct QCD calculations are very compli- 
cated. The model has been used in various approaches 
of the QCD phase transition, both in and out of equilib- 
rium. It has been especially popular in studies involving 
the formation of disoriented chiral condensates. The for- 
mation of DCC was proposed long ago, as a signal of the 
chiral phase transition 40] . 

The C JT method and the linear sigma model have been 
used in other recent investigations, too. Our approach to 
the chiral phase transition appears to be complementary 
in most of these approaches. Throughout this report, 
we have not discussed the issue of renormalization. It 
is well known that the linear sigma model is rcnormaliz- 
able at zero temperature. However, recent attempts to 
incorporate renormalization at finite temperature, have 
not always been successful. An attempt by Amelino- 
Camelia [4^ . to examine renormalization of the linear 
sigma model at finite temperature results in violation of 
the Goldstone theorem. We have mentioned this already 
in the third section. A recent successful treatment of 
the renormalization matters, is the one by Rischke and 
Lenaghan [sij . 

In the Hartree approximation, we have calculated the 
effective masses of sigma and the pions. We call them 
effective masses, since they correspond to the modifica- 
tion of the bare particle masses, which result after their 
"in medium " interactions. Of course, they do not corre- 
spond to the real observed masses, since in our approach, 
we have ignored quantum fluctuations throughout this 
work. 

In the work of Chiku and Hatsuda |5^ |5^ , the opti- 
mal perturbation theory is employed in order to perform 
the calculations. They obtained gap equations for the 
effective masses, which go further than our Hartree ap- 
proximation. They include quantum fluctuations as well. 
So we cannot really compare our results in full. However, 
they do not calculate the effective potential. 

In the work of Nemoto, Naito and Oka [t^, there are 
many similarities to our approach. First of all, the CJT 
method is used for the calculations, however they renor- 
malize the model, and use a definition for the particle 
masses as the minimum of the potential. 

In the work of Caldas, Mota and Nemes |73, the main 
difference is that they do include fermions as well. Their 
calculation of evolution of the condensate look very simi- 
lar to our result. The show a first-order phase transition 
in the chiral limit, and a crossover when there is an ex- 
plicit breaking of chiral symmetry. In the work of Bilic 
and Nicolic |5J|, they also include fermions, and find a 
first-order phase transition. 

Bochkarev and Kapusta have focused in the dif- 
ferences between the linear and non-linear sigma model. 
They use the 0{N) version, and their result is similar to 
our large-iV approximation, since they found a second- 



order phase transition as well. They also find that the 
non-linear sigma model has a second-order phase tran- 
sition as well. 



VI. SUMMARY OF THE RESULTS 
A. Conclusions 

We have studied the chiral phase transition using the 
linear sigma model. In order to understand the nature 
of the phase transition and how it could proceed, in the 
first two sections of this work, we have calculated the fi- 
nite temperature effective potential of this model in the 
Hartree and large- approximations using the CJT for- 
malism of composite operators. This method has the 
advantage that we actually only need to calculate one 
type of diagram. 

In both cases, we solved the system of gap equations 
numerically, and found the evolution with temperature 
of the thermal effective masses. In the Hartree approx- 
imation, we find a first-order phase transition but, in 
contrast, the large A'^ approximation predicts a second- 
order phase transition. This last observation seems to 
be in agreement with different approaches to the chiral 
phase transition based on the argument that the linear 
sigma model belongs in the same universality class as 
other models, which are known to exhibit second order 
phase transitions |4^ . 

However, in the large-A' approximation, the sigma 
contribution is ignored and this, of course, introduces 
errors when we calculate the critical temperature. We 
found that the large- A^ approximation predicts a higher 
transition temperature than the Hartree one. However, 
we have concluded that this is an artefact of the calcu- 
lation, since both Hartree and large-A^ approximations 
should have the same high temperature limit. We should 
recall that, we can calculate the transition temperature 
by considering the high temperature limit from the mass 
gap equation, since it is defined as the temperature where 
the particles become massless. As we have pointed out 
in Section^ in both cases and at high temperatures the 
behaviour of the pion-sigma system approches that of 
the ideal gas. 

When we include the symmetry breaking term ea 
which generates the observed pion masses, we found 
that there is no longer any phase transition, both in 
Hartree and large-A^ approximations. Instead, we ob- 
serve a crossover phenomenon, where the change of the 
order parameter in the Hartree case occurs more rapidly 
in contrast to the smoother behaviour exhibited in the 
large-A^ approximation. Our observation confirms re- 
sults reported recently by Chiku and Hatsuda £^ 
using the ideas of optimal perturbation theory. In their 
analysis they also report indication of a first order phase 
transition in the chiral limit. This observation is closer 
to the real world, since the pions are only approximately 
Goldstone bosons. 
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Of course, as we have already pointed out, the Hnear 
sigma model is only an approximation to the real problem 
which is QCD, but the study of the chiral phase transition 
in the framework of this model could be a helpful guide 
to how one could tackle the original problem, and get 
some insight in the physics involved. For the calculations 
both in Hartree and large-iV approximations, we have 
used the imaginary time formalism which is adequate for 
studies at thermal equilibrium but if one is interested 
in studies of the dynamics of the phase transition, the 
real time formalism seems to be more convenient Q • We 
did not study the system far from equilibrium, however 
the real time formalism would allow us to extend the 
investigation, and consider the dynamics of the system. 
In real time formalism, we have the advantage that the 
propagator splits into two parts from the beginning, so 
making it easier to calculate thermal corrections to the 
effective masses. 

However, our attempt to go futher than the Hartree 
approximation has run into serious difhculties. It is not 
absolutely clear which class of diagrams we should con- 
sider. Our motivation to include only the sunset-type 
diagrams with two thermal particle only, comes out from 
conciderations of the low-energy theorem, and our per- 
cistence into symmetry principles. This is why we have 
not tried to renormalise the model as well. It is well 
known that the linear sigma model is a renormalizable 
theory at zero temperature, and that finite temperature 
effects do not indroduce new divergencies - at least at 
the ultraviolet -, since the temperature acts as a natural 
cut-off. However, renormalization of the model is inves- 
tigated in other recent investigations as it was mentioned 
already. 

The fact that near the transition temperature, the ef- 
fect of higher loops may become important, it is well 
known for a long time, and suggests that further investi- 
gation is needed. We would suggest that an extension of 
our two loop calculation is tractable, provided one takes 
into account the efffects of all the loops contributing at 
the two-loop level, as a first step. However, it is not clear 
how the C JT method could be used to include the effects 
of higher loops. Selfconsistent approximations have been 
used for a long time, but it is not clear how they work 
beyond the one loop level. 



B. Some recent results 

We conclude this work reviewing some recent progress 
on the subject. Since the work which presented in the 
previous sections was initially written, there has has been 
quite some progress on similar studies. The linear sigma 
model has always been a very popular model in studies 
attempting to mimic the most of the low energy region 
of QCD. 

Perhaps the most interesting result, which is related 
to this work and has been published recently, is the work 
of Baacke and Michalski ||8J] , where the authors examine 



the 0{N) linear sigma model beyond the Hartree ap- 
proximation. The same authors have also studied the 
sigma model at non-equilibrium in (see also [s^'). 
However, it seems that their approach does not suffer 
from the difficulties that we have faced trying to solve 
the full system of gap equations at the two loop level of 
the effective potential. An essential difference between 
our approach and the one by Baacke and Michalski is the 
renormalisation of the model. As we have seen already 
in Section the CJT method is based on a resumma- 
tion of the 2PI graphs. The analysis of [s^l is based on 
the a variation of the CJT method, which is called 2PPI 
resummation. This method has been initiated by Ver- 
schelde and C opp ens [87l | and it is presented in a number 
of papers [8l|8g,'9S]. 

In order to study the effective potential Baacke and 
Michalski have calculated the full subset of sunset dia- 
grams. This of course makes their approach more com- 
plete. However, our attempt was focused mainly on the 
sunset with two thermal pions trying to test if the low en- 
ergy theorem has an analogy at finite temperature, rather 
than tackling issues like for example renormalization of 
the model at the full two loop calculation, or attempts 
to consider the full two loop case. 

An alternative approach to the problem of the selfcon- 
sistent approximations, is the work of Knoll and van Hees 
[Olt .92- .93,) where they use the $ derivable approxima- 
tion. This method orinates from the work of Baym 
and it has essentially many similarities with CJT. How- 
ever, it seems that renormalisation is controlled better 
under this calculational scheme. 

The CJT method is an attractive framework in order 
one to perform selfconsistent approximations and has at- 
tracted quite some attention. The method is used in a 
recent paper by Roder et al [o^l where the authors study 
restoration of chiral symmetry and the low lying meson 
spectrum at finite temperature. Also in another recent 
study of the 0(4) linear sigma model Phat et al 
have used the CJT method as well. In their analysis, 
in contrast to our approach, they have studied the issue 
of renormalisation while they have incorporated the ef- 
fects od higher loops. Other recent work on the 0{N) 
linear sigma model, but without using CJT, is the work 
of Patkos et al J)6, 97, 98, 99] where they study the be- 
haviour of the spectal functions, as well as, the position 
of the sigma pole on the complex energy plane. 
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The result is summarised as 
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APPENDIX B: ONE LOOP MASS CORRECTION 



APPENDIX A: THE POTENTIAL AT ONE-LOOP 

In this appendix, we review the calculation of the one- 
loop contribution to the effective potential for a At/)^ the- 
ory, using the imaginary time formalism. This is a well 
presented calculation, and appears in many research pa- 
pers and textbooks, as for example references [TMl20ll2^ . 
However, in order to make this presentation more self- 
contained, and also compare it with the real time cal- 
culation, we will reproduce the basic steps here. The 
propagator is given by 



(Al) 



The one loop contribution to the effective potential is of 
the form 



ln(fc'^ -I- m^) 



(A2) 
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where E"^ = k^ + m^, and we have indicated the sum 
over the Matsubara frequency fco = ujn — 27rn/3^^ ex- 
plicitly. We should recall that f3 = 1/T, where T is the 
temperature. 

In order to evaluate this sum we can define 
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so the derivative is 
du{E) 

dE " ^ 4^2^2//3^ + E"^ 

n ' ' 

We can now use the formula ^ 
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We have seen in n Section that the one-loop correc- 
tion to scalar propagator only causes a mass shift. The 
relevant graph is shown in Fig. 1331 . 




FIG. 33: The self-energy contribution to scalar propagator. 

This graph involves evaluation of the following integral 
d^k 1 

(Bl) 
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If we use the formula 
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we find that 
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From the fact that 

coth/3a; = 1 + 2n(2x) , (B4) 
where n{E) is the Bose-Einstein distribution function 
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e/3£ _ 1 ' 

we find that he above integral is written as a sum of two 
parts 

F(m) = Foim) + Fp{m) (B6) 
1 f d^k I 1 r d^k I 1 



2 J (27r)3 2E 2j (27r)3 E ef^^ - 1 
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APPENDIX C: THE KELDYSH CONTOUR 

The imaginary time or Matsubara formalism has the 
disadvantage that if one wants to calculate non equilib- 
rium quantities, an analytic continuation to real time 
must be performed. In this formalism, we are integrat- 
ing in the complex plane of the Minkowski time, from 
t = to t — ip. The periodicity of the fields enables us to 
generalise this interval to, Iq to to — if], for any real to. 

Alternatively, one can work directly in real time. In 
this case, we need a different choice of contour with the 
same endpoints which will contain the real time axis. A 
choice of contour for time integration, which is called the 
Keldysh contour, is shown in Fig. |^ 
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FIG. 34: Keldysh contour 

As we can see in Fig. |M1 the Keldysh contour consists 
of four pieces. However, there is an alternative choice 
of contour [^. For analyticity reasons, it can be shown 
that the contributions from C3 and C4 can be neglected. 
Therefore, we arc left with two possibilities for both the 
initial and final points of integration to lie either on Ci 
or C2. This gives four different propagators, which are 
usually written in matrix form. We have given these 
propagators in Section llDl 



APPENDIX D: 



THE PION SELF-ENERGY 
GRAPHS 



1. The pion self— energy graph with one thermal 
pion 

Although the calculation of the self-energy graphs for 
A(/)^ is a well known subject, and appears in many re- 
search papers, we repeat the derivation here, hoping to 
make the presentation of the calculations more complete. 
Evaluation of the self-energy can be done using any of 
the formalisms: ir nag inar y t ime, real time or thermo- 
field dynamics |2ll ISOl ISU 183 . l83l | . We find the real time 




FIG. 35; The pion self-energy diagram with one thermal pion. 



convenient at this stage. The pion self-energy graph with 
one thermal pion in the loop, is given in Fig. 1351 

As mentioned already in Section IIVI in our case, we 
need only the (1,1) component of the matrix propagator. 
This has a zero temperature part 



— m + It fc — 

(Dl) 

where P means the principal value, and a thermal part 
given by 



Vrik, T) = 27rJ(fc2 - rT?)n{ko) 



(D2) 



where niko) is the Bose-Einstein distribution function 

1 



n{ko) 



exp(/3|fco|)-l ' 
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with (3 — 1/T the inverse temperature and fcg = 
•\/k^ + irfi . 

We use the following convention for 4-momenta = 
(fco,k), and reserve k for the magnitude of the 3- 
momentum, so fc = |k|. Since we are interested in the 
real part of the self-energy graph, we ignore the second 
part of the vacuum propagator (the one with the delta 
function). Then the exact expression corresponding to 
the graph in Fig. 1351 is written as 
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As we have mentioned earlier, the superscripts denote to 
which particle's self-energy we are reffered to, while the 
subscripts denote which particles are running in the loop. 
A subscript 9 denotes a thermal propagator. 
By using partial fractioning, we can write the above in- 
tegrand as 
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Where the superscript in the expression for the energy 
has the obvious meaning, E^ = yj + m^. In order to 
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perform the integration, we use the following property of 
the Dirac delta function 
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Inserting everything in the integral, and integrating over 
the delta function, results in 
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or, in a more compact form, as 
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In the case when we consider pions as static, then p = 
and so po = M^^ , therefore we can write the above as 
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or even as 
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2. The pion self-energy graph with one thermal 
sigma 
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FIG. 36: The pion self-energy diagram with one thermal 
sigma. 



now, this graph is absolutely analogous to the one pre- 
sented above. Therefore, following the same procedure 
as before, we end up with the following expression 
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APPENDIX E: THE SIGMA SELF-ENERGY 
GRAPHS 

1. The sigma self— energy graph with one thermal 
sigma 

Evaluation of the sigma self-energy graphs can be done 
in a completely analogous way as to one of the pions. 
However, since it involves identical particles it is even 
easier. Therefore, 
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FIG. 37: The sigma self-energy diagram with one thermal 
sigma. 



Another graph which contributes to the pion self- 
energy, is the one given in Fig. 1361 As we can easily 
observe, apart from that the thermal particle is a sigma 
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2. The sigma self— energy graph with one thermal 
pion 

In this case we have the topology as in Fig. OHl 
Evaluating the relevant integral as above, we find 




FIG. 38: The sigma self-energy diagram with one thermal 
pion. 
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APPENDIX F: CALCULATION OF SUNSET 
GRAPHS 

1. The sunset with two thermal pions 

The sunset diagrams consist of a set of diagrams with 
the topology given in Fig. 1241 . However, in the approx- 
imation adopted for this work, we ignore quantum fiuc- 
tuations, so we do not take into account the graphs with 
one thermal propagator or none. The graphs with three 
thermal propagators do not contribute (because of con- 
servation of energy), so we are left only with the graphs 
containing two thermal propagators. 




FIG. 39: The sunset diagram, with two thermal pions. A 
line with a dash corresponds to the thermal part of the full 
propagator. 

In the case of two thermal pions in the sunset as shown 
in Fig. 123 we deal with the expession 

2Trn(ko) ztt n[po) 



(27r)4 J (27r)4- 



M2) 



(P + Kf - Ml 



(Fl) 



and in order to evaluate this integral, we use the prop- 
erty of the delta function given by Eq. Then our 
expression will consist of four terms containing products 
of delta functions of the form 



5{ko±E^)5(po±E^) 



(F2) 



In order to perform this integration, it is convenient to 
choose k as defining the polar axis, so then 9 is the angle 
between the three vectors k and p. Then, the first term 
of the integral given in Eq. HF1|I will appear as 



RaiM^,M,) 



1 



dko kdk 



2(27r)4 J ko €0^0 - 1 
dpo pdp 



where we have defined A — A{k,p, 9) as a shorthand of 
the lengthy expression 

kpsme d9 5{ka - E^) S(pn - Ep) 



A 



kl~ B +pI- p"^ + 2kQpo - 2kp sin 9 - M^ 



After performing the delta function integration, we end 
up with the result 



Ra(M^,M, 



1 



kdk 1 



2(27r)4 J E^ el^^i - 1 
pdp 1 



El e^^? - 1 



kpB{k,p,9) 



and we have defined again a shorthand expression B 
B{k,p, 9) as 



B 



kp sin 9 d9 



{E'^ + ElY 



2fcpsine'- Af2 



Finally, performing the angular intergration, wc find 



RaiM^,M„) = 



2(2 



X In 



ln)^J E^ ^ J El ^ ^' 



2 Ml ~ Ml + 2ko po + 2kp 



2 Ml- Ml + 2ko Po ~ '2kp ' 

In repeating this procedure, there will be four other simi- 
lar terms, one for each of the combinations of delta func- 
tions given by Eq. ljF2p . Then the final result can be 
written as 



{M,,M,) 



1 



kdk 
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where 
A{k,p,M^,M^) ^ 

B{k,p,M^,M,) = 



Ml{e+p^) + Ml(Ml 
+ {2Ml - Ml)kp , 

Ml{e + p^) + Ml{Ml 
-{2Ml-Ml)kp. 



(F3) 
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2. The sunset with two thermal sigmas 

In the case when we need to calculate the sunset with 
two thermal sigmas, the situation will be absolutely anal- 
ogous to the one presented above for the two pions. The 
diagram has the topology given in Fig. 021 ■ 



■ p 

P + K 



a 
K 



FIG. 40: The sunset diagram, with two thermal sigmas. A 
line with a dash corresponds to the thermal part of the full 
propagator. 



Then the exact expression corresponding to this graph 
is given by 



27r)4 J (27r)4 
S{K^ - Ml) S{P^ - Ml 



27m(fco) 27rri(po) 
(F4) 



[P + KY ~ Ml 

In order to evaluate this, we repeat the steps as in case 
with the two thermal pions. The final result will be of 
the form 



laaeae {M^,M^) 



: In 



4(fc+p)2 + 3M2 + 4fcp 
4(A;+p)2-|-3M2~4fcp 



(F5) 
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